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Abstract. We present Backlund transformations (BTs) with parameter for 
certain classical integrable n-body systems, namely the many-body generalised 
Henon-Heiles, Garnier and Neumann systems. Our construction makes use of 
the fact that all these systems may be obtained as particular reductions (sta- 
tionary or restricted flows) of the KdV hierarchy; alternatively they may be 
considered as examples of the reduced sl{2) Gaudin magnet. The BTs pro- 
vide exact time-discretizations of the original (continuous) systems, preserving 
the Lax matrix and hence all integrals of motion, and satisfy the spectrality 
property with respect to the Backlund parameter. 



Backlund transformations (BTs) are an important aspect of the theory of inte- 
grable systems which have traditionally been studied in the context of evolution 
equations. However, more recently there has been much interest in discrete sys- 
tems or integrable mappings ||l|, |l^. Within the modern approach to separation 
of variables (reviewed by Sklyanin in flql) this has led to the study of BTs for 
finite-dimensional Hamiltonian systems [|l2[ . The latter are canonical transforma- 
tions including a Backlund parameter A, and apart from being interesting integrable 
mappings in their own right they also lead to separation of variables when n such 
mappings are applied to an integrable system with n degrees of freedom. The se- 
quence of Backlund parameters Xj together with a set of conjugate variables iij 
constitute the separation variables, and satisfy a new property called spectrality 
introduced in 

We proceed to develop these ideas with some new examples of BTs for n-body 
systems, namely the many-body generalisation of the case (ii) integrable Henon- 
Heiles system, the Garnier system and the Neumann system on the sphere (see Q ) . 
It is known that the case (ii) Henon-Heiles system is equivalent to the stationary 
flow of fifth-order KdV ||] , while the Garnier and Neumann systems may be ob- 
tained as restricted flows of the KdV hierarchy ||T^ . Thus we derive BTs for these 
systems by reduction of the standard BT for KdV, which arises from the Darboux- 
Crum transformation |^ for Schrodinger operators. The restriction of the Darboux 
transformation to the stationary flows of the modifled (mKdV) hierarchy has been 
discussed in j^. 

In the following section we describe how the reduction works in general, before 
specialising these considerations to each particular system and presenting the asso- 
ciated generating function for the BT. We note that these systems are examples of 
the reduced Gaudin magnet , so that we have the following Lax matrix 



where (up to scaling) the Sj satisfy n independent copies of the standard sl{2) 
algebra: 



1. Introduction 



(1.1) 




(1.2) 



{SlS^} = ±26,kS, 



{S 
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For the Hcnon-Heiles and Garnier systems the matrix B(u) is respectively qua- 
dratic and hnear in the spectral parameter u, while for the Neumann system it 
is independent of u and turns out to be constant due to the constraint that the 



particles lie on the sphere (hence the Poisson algebra (1.2) must be modified by 
Dirac reduction). 

We have constructed the BT for the (non-reduced) sl{2) Gaudin magnet with 
quasi-periodic boundary condition in |l^, while some preliminary results on the 
classical Garnier system and two-body Henon-Heiles system first appeared in 

2. Classical integrable systems and KdV 

2.1. Restricting the BT. As is well known, the Darboux-Crum transformation 
m consists of mapping the Schrodinger operator df + V — X to another operator 
df + V — Xhy factorizing the former and then reversing the order of factorisation. 
Given an eigenfunction (jj satisfying 

{d^ + F - A)0 = 

we may set y — (log[(/)])t and then 

(2.1) V = ^yt^y^ + X, V ^ y^ - y^ + X; 

for A = this is just the Miura map for KdV. Also given another eigenfunction ^ 
of the Schrodinger operator with potential V for a different spectral parameter u 
we have 

where the transformation to the new eigenfunction and its derivative may be 
given in matrix form as 



for any constant k. From (^J|) follows the standard formula for the Darboux- 
Backlund transformation of KdV, V = V + 2{log[(j)])tt- 

For what follows it will also be necessary to consider a product of eigenfunctions 
for a Schrodinger operator with potential V and eigenvalue u, 

with Wronskian "0*^' — V'V't — 2771. It is well known that / satisfies the Ermakov- 
Pinney equation 

(2.3) fftt - Iff + 2{V - u)f + 2m2 = 0. 



If we now transform -0 and ■0' according to ( |2.2| ) then we find a new product of 
eigenfunctions / = "00' satisfying the same Ermakov-Pinney equation but with V 
replaced by V , given explicitly by 

(2.4) /^(A-,,)-ii^!^, Z=\u^yf, 

where we have set — {X ~ u)^^ to ensure that the transformed eigenfunctions 
have the same Wronskian 2m. It is also straightforward to show that, in terms of 
/, the quantity Z can be written as Z = — ^/t — yf (see |^). 

We can now describe how this transformation restricts to the finite-dimensional 
Hamiltonian systems presented below. The systems are expressed in variables 
{qj,Pj) which appear in the Lax matrix (LI) via the identification [pT| Q 
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Sj =Pjqj, Sj ~ ^Vj ^ 2"' ~ '^j ■ 

9i" 
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For Hcnon-Heiles and Garnicr the non- vanishing Poisson brackets are the standard 
ones {pj, Qk} = Sjk which provide a reaUzation of the algebra (1.2); for the Neumann 



system on the sphere the bracket must be modified by Dirac reduction. 

All of the systems are Liouville integrable, and thus have a complete set of Hamil- 
tonians in involution, but for these purposes we concentrate on the Hamiltonian 
h generating the flow corresponding to t above (in KdV theory this is usually de- 
noted X, the spatial variable). For this flow the Lax equation Lt = [N,L] is the 
compatibility condition for the linear system 

(2.5) i(u)* = t;*, ^t^N^-, ^""(l """^J*'^^'^ 

Observe that the second part of the linear system is just a Schrodinger equation 
for the potential V; for Neumann and Garnier this is a function of {qj,pj) for 
J — 1, . . . , n, while for Henon-Heiles there is an extra pair of conjugate variables 
(qn+ijPn+i) such that V = q„+i. 

The equations of motion generated by this Hamiltonian take the form qj^t — Pj 
and 

2 

TTl - 

(2-6) Pj,t = qj,tt = (flj - V{qk,Pk))qj - -4 

"j 

for j = 1, . . . , n; for Henon-Heiles there are also equations for q„+i and Pn+i = 



qn+i.t- The important thing to observe is that (2.6) is equivalent to the fact that 
— q'j satisfies the Ermakov-Pinney equation (2.3) corresponding to a Schro- 
dinger equation with potential V and eigenvalue aj. Thus to obtain a Backlund 
transformation for these many-body systems we simply apply a Darboux-Crum 
transformation ( |2.lD to the potential V = V{qj,pj) to obtain V = V{qj,pj), and 
then we know that the solutions of the Ermakov-Pinney equation must transform 
according to ( |2.4| ). By this procedure we may explicitly construct the BT for the 
many-body systems below (or for any restricted flow of KdV) , and it is then simple 
to calculate the generating function F[qj, qj) of this canonical transformation, such 
that 

dF = ^{jpjdq.j -pjdqj). 

3 

The discrete Lax equation for the BT, 

ML = LM 

where L — L{qj,pj;u), is necessary to ensure the preservation of the spectral curve 
det(u — L{u)) = (so that all the Hamiltonians in involution are preserved). This 
follows immediately from the properties of the Darboux-Crum transformation, since 



we know that the vector ^ in the linear system (2.5) must transform according to 



(2.2), and hence we may take (setting k = I) 

y y'^ + u~ X 

-y 



(2.7) M - ^ 



Of course we must determine y as a function of the dynamical variables. In the Gar- 
nier and Henon-Heiles cases it turns out that the potential depends on coordinates 



only, V — V{qj), and so by adding the two equations in (2.1) we obtain 



yiq,,q~,) = ±^X-^iV + n 

to obtain the correct continuum limit of the discrete dynamics it is necessary to 
take the negative branch of the square root (see |§, |§]). For the Neumann system 
V depends on both coordinates and momenta, so the above does not yield y{qj,qj). 
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There is another way of writing L which arises more naturaUy via reduction from 
the zero curvature representation of the KdV hierarchy ^ viz 



L{u) - ' ^ ' 2 



n -in* 



where 

(2.8) n(^) = 5]^^ + A(z.). 



9? 



A is a polynomial in u fixing the dynamical term B in (f.l); we shall present the 



appropriate A and B in each case below. Clearly the t derivatives of 11 can be 



rewritten using the equations of motion to yield (1.1) 



Finally if we write the (hyper-elliptic) spectral curve as 

= R{u) 

then it is easy to check that the spectrality property is satisfied for these 
systems, in the sense that defining the conjugate variable to A by 

dF 

we find that 

with eigenvector = (y, 1)^, or in other words [i} — R{X) so that (A,/i) is a point 
on the spectral curve. Note that (as for the examples in [^0|, |l2|) this eigenvector 
spans the kernel of M, 

M{X)n = 0. 

We can also write y explicitly in terms of both the old and the new variables related 
by the BT, thus: 

, , n.(A) + 2A. (nt(A)-2A.) 

(2.9) yiq^,p^) = —-^, y{q,.P,)- ^5^; 

clearly we denote n(A) = W[qj,p.j \ A). 

2.2. Generalised Henon-Heiles system. For the many-body generalisation of 
case (ii) integrable Henon-Heiles system, the Hamiltonian generating the t flow 
takes the form 

1 /l " \ 1 " / 

The original case (ii) integrable Henon-Heiles system corresponds to n = 1 with 
c = rrij — Qj = 0. The link between stationary fifth-order KdV and the type (ii) 
system was noted by Fordy in |^ , although this was anticipated in work of Weiss 
|l8| , who used Painleve analysis to derive a BT and associated linear problem (a 
similar result also appears in ||l^ ) . None of these authors wrote a BT explicitly as 
a canonical transformation with parameter, although (without parameter) this was 
done for a non-autonomous version in . 

For the Lax matrix L of the generalised (n + l)-body Henon-Heiles system we 
have A = —16u — 8qn+i so that the extra term -8(11) is given by 

^ = ( -mu'-sL, 4p1, ) ' ^ = -16"^ + - 4.?.,r - E .| " 4c. 
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The equations of motion for h imply that the squares of the first n coordinates g| 
satisfy the Ermakov-Pinney equation ( |2.3| ) for m — rrij with 

V = Qn+l 

and eigenvalue aj. Thus the BT for the system can be calculated directly by 
applying the Darboux-Crum transformation to = Qn+i, to yield V = (jn+i, and 
applying (2.4) to each ior j ~ 1, . . . , n. 

After some calculation the generating function for this canonical transformation 
is found to be 

16 



n9,.%;A) = 5:(z, + ^iog 
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K+i + 2g„+ig„+i + 2g2+i + i J2{q^ + «i ) + 2cj y, 

where we have found it convenient to use the quantities Zj(qj,qj) and y{qj,qj) 
defined by 

(2.10) Z] ^m^^ + iX-a,)q]ql 

and 

y = - ^{qn+i + qn+i)- 

In order to check the spectrality property, we have explicitly found that the 
eigenvalue of L{\) with eigenvector fl = (y, 1)^ can be written as 

, ~ u >^ Z, IdF 

^ A - y dy 

which precisely equals —2^ as required. 

2.3. Garnier system. For the Garnier system the t flow is generated by the Hamil- 
tonian 



] = 1 ] = l 3 = 1 \ ^3 I 



This differs from the traditional Garnier system as in ^, |l^, |l9| by the inclusion of 
extra inverse square terms. The Newton equations for the qj are 

^2 



^Z^Qkh] = ajqj - -3- 



so clearly for the standard restricted flows of KdV |]T9j, when mj = 0, each qj is an 
eigenfunction of a Schrodinger operator with potential 

^ = 2^,1 

j 

and eigenvalue aj, while in general is a product of eigenfunctions satisfying the 
Ermakov-Pinney equation Q for m = mj . 

The Lax matrix of the Garnier system has A = 1, so i takes the form (1.1) with 

^^(0 u-Z^q^ 



1 

Applying the Darboux-Crum transformation we obtain a new potential 

3 
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and the corresponding BT induced on the Garnier system is equivalent to gauging 



L by the matrix M of the form (2.7) with 



Finally we can calculate the generating function for this BT, which may be 
written as follows: 



F{qj,qj;X) 



■log 



1 



where y{qj,qj) is as above and Zj is given by the same expression (2.1C) as for 
Henon-Heiles. In ||] we derived this generating function for the special case nij = 
when the logarithm terms do not appear. To check spectrality we notice that L(X) 
has eigenvalue 



n{qj,qj;X) 



^ A 



with eigenvector f2, and so we see that /i 



2.4. Neumann system on the sphere. For the Neumann system the t flow is 
generated by 

i2 ' 




Once again this has extra inverse square terms compared with the standard Neu- 
mann system jl^ , . The Poisson bracket for this system is modified by constrain- 
ing the particles to lie on a sphere, so that 



(2.11) 



(9,9) 



g| — const, 



3 



qjPj = 



which results in the non-vanishing Dirac brackets 



(2.12) 



{pj,qk} = Sjk - 



QjQk 



{PJ^Pk} 



QjPk - QkPj 



iq,q)' ^"^'"^^ iq,q) 
With this bracket the Hamilton equations are qj^t — Pj and 

V=iq,qr^Y.(p3+-.^^-4]- 



with 



The Lax matrix for the Neumann system arises by setting A = 0, which in (1.1) 
gives the following matrix B: 

{q,q) 




B = 



In fact if we start from the linear system (2^) and leave V unspecified then (2.6) as 
well as the constraint {q, q)t = are consequences of the Lax equation, and together 
these are sufficient to determine the form of V\ this is also how the equations for 
the constrained Neumann system arise in a Lagrangian approach jl^. 



Given that the phase space is now degenerate with two Casimirs given by (2.11), 
it would appear that the standard sort of generating function will no longer be 
appropriate for describing a BT. It turns out that we can apply the Darboux-Crum 



transformation as before, and transform the quantities q^ according to (2.4). In this 



way we obtain new variables qj{qk,Pk) and Pj{qk,Pk), which are natu rally written 
with the use of the quantity y{qk,Pk) given by the first formu la i n (2.9); on the 
Lax matrix this transformation arises by gauging with M as in (2.7). Similarly the 
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transformation can be inverted to give qj{qk,Pk) and Pj{qk,Pk) written in terms of 
yiqk,Pk) given by the right hand formula of ( |2.9| ). 

However, it would still be nice to write a generating function for this transfor- 
mation. We have found that if we formally take 



2 



F(q„g,;A)=5] Z, + -^log 
with Zj given by ( 2.10| ) as usual, and regard y as a sort of Lagrange multiplier 



(independent of the coordinates and A), then we do indeed obtain the correct ex- 
pressions 

_ dF _ dF 

- - ~w 

but these contain y which is unspecified. If we then require that the constraints 



(2.11) are preserved under the BT applied from old to new variables or vice- versa, 
then in either direction the constraints are preserved if and only if y satisfies a 
quadratic equation with solution given respectively by the formulae (|2.9| ). Alterna- 
tively if we require spectrality then second component of the equation L{X)n — /iil 
gives 

. ~ ^^ \^ Z, dF 

.7 = 1 



as required, while the first component gives (after making use of the formula (2.10) 
and the BT) 

Hence spectrality requires that the second term vanishes, and so the first constraint 
( ^.11 ) is preserved; the preservation of the second constraint is then an algebraic 



consequence of the BT. 

Thus we see that for this BT we can write the new variables as functions of the 
old and vice- versa, but a formula for y{qj,qj; X) is lacking. Also this discretization 
of the Neumann system is apparently new, since it is exact (preserving the Lax 
matrix for the continuous system) unlike the Veselov or Ragnisco discretizations 
discussed in p4| . 

3. Conclusions 

It would also be interesting to look at BTs with parameter in the non-autonomous 
case [0 , where deformation with respect to the Backlund parameter would probably 
have to be introduced (corresponding to the associated isomonodromy problem). 
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